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ABSTRACT 

Knowledge of the mass distribution within galaxies is essential for our understanding of their structure, but is 
challenging to acquire due to the complex spatial distribution of stellar populations. The advent of deep, high- 
resolution data over a wide range of wavelengths provides an unprecedented opportunity to measure spatial 
variations in stellar populations and quantify their effects on the mass-to-light distribution. In this Paper we 
present stellar mass surface density profiles of a mass- selected sample of galaxies at < z < 2.5, obtained using 
very deep HST optical and near-infrared data. This is the first time accurate stellar mass surface density profiles 
have been measured for galaxies up to z = 2.5. The profiles are constructed using radial rest-frame u-g color 
profiles, from which we derive the stellar mass-to-light ratio as a function of radius. The key advantage of this 
study compared to previous work is that the profiles are deconvolved for PSF smoothing, allowing accurate 
measurements of high-redshift galaxy structure. The profiles are obtained using a measurement method which 
enables us to trace complex galaxy structures such as rings and faint outer disks. We derive stellar half-mass 
radii and find that these are on average ~ 25% smaller than rest-frame g band half-light radii. The average 
size difference is the same at all redshifts, and does not correlate with stellar mass, specific star formation rate, 
effective surface density, Sersic index, or galaxy size. Although on average the difference between half-mass 
size and half-light size is modest, for approximately 10% of massive galaxies this difference is more than a 
factor two. These extreme galaxies are mostly extended, disk-like systems with large central bulges. 

Subject headings: cosmology: observations — galaxies: evolution — galaxies: formation — galaxies: high- 
redshift 



1. INTRODUCTION 

Over the past decades quantitative studies of high-redshift 
galaxy structure have advanced tremendously. Sensitive, high 
resolution instruments such as the Hubble Space Telescope's 
(HST) Advanced Camera for Surveys (ACS) and the Wide 
Field Camera 3 (WFC3) have made it possible to measure the 
structure of faint high-redshift galaxies at sub-kpc scales. Fur- 
thermore, the availability of easy-to-use photometric redshift 
and stellar population fitting packages has made it possible 
to straightforwardly measure a wide variety of parameters for 
ever-increasing numbers of galaxies. 

Since a small amount of recent star formation can have a 
disproportionately large contribution to a galaxy's light com- 
pared to its mass, galaxies are usually observed at the redmost 
wavelengths, where emission from young stars is weakest. At 
low redshift, this can be done quite effectively, since rest- 
frame near-infrared (NIR) data is available at high enough 
resolution. At higher redshifts, however, it is impossible to 
observe at such long wavelengths with sufficiently high an- 
gular resolution. Until recently, the Hubble Space Telescope 
(HST) Advanced Camera for Surveys (ACS) was the only 
wide-field instrument capable of measuring the structure of 
high-redshift galaxies in any detail. At z = 2, the reddest fil- 
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ter available on ACS, Zgso, corresponds to rest- frame near- 
ultraviolet (NUV) wavelengths. Use of such short- wavelength 
data has been shown to result in drastically different conclu- 
sions about galaxy structure and morphology, compared to 
rest- frame optical data (e.g ., Labbe et"aL||2003[ Toft et al. 
[20051 |Cameron et al.|201 

With the introduction of the HST Wide Field Camera 3 
(WFC3) it has become possible to measure rest-frame optical 
light of z ~ 2 galaxies at a resolution approaching that of the 
ACS. The redder light detected by this instrument provides a 
much better proxy for stellar mass. However, color gradients 
are known to exist to some extent in all types of galaxies at 
redshifts up to at least z ~ 3, such that most gal axies contain 
a relatively r ed core and blue outer regions (e.g., [van Dokkum 
[eTaLl|20T0l |Szomoru et aLl|20TTl |Guo et al.||Z01i| ). These 
color variations are caused by a combination of varying dust 
content, metallicity and stellar age, and imply that the stellar 
mass-to-light ratio (M*/L) of a given galaxy varies with po- 
sition within that galaxy. Thus, even though rest-frame opti- 
cal light is a better tracer of stellar mass than rest-frame NUV 
light, neither accounts for the complexity of stellar population 
variations within galaxies. 

By fitting stellar population models to resolved galaxy pho- 
tometry, it is in principle possible to infer spatial variations 
in stellar mass, age, metallicity, dust content, and other pa- 
rameters. This approach is currently somewhat limited by the 
lack of high-resolution data at infrared wavelengths, but can 
nonetheless be used to measure several basic stellar popula- 
tion properties. An example of this technique is presented in 
Wuyts et al.| (2012), who have performed stellar population 



modeling on resolved HST data, using integrated IR observa- 
tions as constraints on the overall properties of their galaxies. 
In this approach the integrated photometry serves as an im- 
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portant tool to constrain the overall spectral energy distribu- 
tion (SED) of a galaxy, while the HST data provide informa- 
tion regarding the spatial variation of the stellar populations 
within these constraints. 

In this Paper we explore an alternative method to recover 
M*/L variations, using a simple empirical relation between 
rest- frame u-g color and M*/L. Using this method we con- 
struct stellar mass surface density profiles corrected for the 
effects of the PSF, for a mass- selected sample of galaxies 
between z = and z = 2.5. We compare the resulting half- 
mass radii to half-light radii based on rest- frame optical imag- 
ing. All sizes presented in this Paper are circularized sizes: 
r e = re,a\Jb/a. Throughout the Paper we assume a ACDM 
cosmology with Vt m = 0.3, = 0.7, and H Q = 70 km s -1 
Mpc- 1 . 

2. DATA AND SAMPLE SELECTION 

2.1. HST imaging 

We make use of deep near-IR imaging of the GOODS- 
South field, obtained with HST/WFC3 as part of the Cos- 
mic Assembly N ear-infrared Deep Extragalactic Legacy Sur - 
vey (C ANDELS ; |Grogin et al.pOTT] [Koekemoer et al.|2QTT| ). 
When completed, this survey will cover ~ 700 arcmin 2 to 2 
orbit depth in 7 814 , J l25 and H m (COSMOS, EGS and UDS 
fields), as well as ~ 120 arcmin 2 to 12 orbit depth (GOODS- 
South and GOODS-North fields). We use the deepest cur- 
rently available data, which consist of 9 orbits in /125 and 
//160 taken over GOODS-South. These NIR data are com- 
bined with deep HST/ ACS data in the #435, V606> ^775 and 
Zg5o filters from the Great Observatories O rigins Deep Survey 
(GOODS ACS v2.0; |Giavalisco et al.|2004l ). The full width at 
half-maximum (FWHM) of the point- spread function (PSF) 
is w 0.12-0.18 arcsec for the WFC3 observations and ^0.1 
arcsec for the ACS observations. The WFC3 and ACS im- 
ages have been drizzled t o pixel scales of 0.06 an d 0.03 arcsec 
pixel -1 , respectively (see Koekemoer et al.|20lT fo r a detailed 
description of the C ANDELS data reduction, and |Giavalisco 
|et al.|2004| for details of the ACS reduction). 

We select galaxies usi ng the ^-selected FIREWORKS cat- 
alog (Wuyts et al. 2008). This catalog combines observa- 
tions of the Chandra Deep Field South (CDFS) ranging from 
ground-based £/-band data to Spitzer 24/im data. It includes 
spectroscopic redshifts where ava ilable, as well as phot omet- 
ric redshifts derived using EAZY (Brammer et al. 2008). The 
photometric redshifts have a median Az/(1 Tz) = -0.001 with 
a normalized medi an absolute deviation of cf^mad = 0.032 
( Wuy ts et al.j[2 008 ). Stellar masses are estimated from SED 
fits to the full photometric data set (N.M. Forster Schreiber 
et al. 2012, in preparation), assuming a Kroupa IMF a nd the 
stellar population models of Bruzual & Chariot ( 2003). Star 
format ion rates have been calculated using the U V and 24 fim 
fluxes ( |Wuyts et al.|20Q9j ). 

We limit our analysis to galaxies with redshifts between 
z = 0.5 and z = 2.5; within this wavelength range we have 
the wavelength coverage needed to measure rest- frame u and 
g band photometry. We select galaxies with stellar masses 
above 1Q 10 - 7 Mq, which is the completeness limit in this 



2.2. SDSS imaging 

We compare our high-redshift galaxies to low-redshift 
galaxie s observed as part o f the Sloan Digital Sky Survey 
(SDSS; Abaza jian et al.|2009| ). In order to obtain the deepest 
possible galaxy photometry we li mit our analysis to data from 
SDSS stripe 82 ( Annis et al.|2011| ). This region of the sky has 
been repeatedly imaged by SDSS, resulting in data which are 
~ 2 magnitudes deeper compared to standard SDSS imaging. 
The PSF FWHM is « 0.6 arcsec, and the images have a pixel 
scale of 0.396 arcsec pixel -1 . 

We base our g alaxy selection on th e NYU Value-Added 
Galaxy Catalog (Blant on et al.|[2005] ); stellar masses and 
star formation rates are ta ken from the MPA-JHU catalog^] 
(Brinchm ann et al.]|2004| ). We select galaxies with spec- 
troscopically measured redshifts in a narrow redshift range 
z = 0.06 ± 0.005 and with stellar masses above 10 ia7 M© . This 
results in a sample of 220 galaxies. 

3. ANALYSIS 

3.1. Rest-frame surface brightness profiles 

Most studies of galaxy structure at high redshift are based 
on parametrized two-dimensional surfac e brightness profil e 
fits (e.g., using the GALFIT package of jPeng et aLl |2002). 
However, such methods do not to account for deviations from 
the assumed model profile, which is generally an r l l n Ser- 
sic profile. The techniq ue used in this paper , which is de- 



redshift range (W uyts et al.||2009| ). This redshift and mass 
cut results in a sample of 183 galaxies, of which 110 are at 
0.5 < z < 1.5 and 73 at 1.5 < z < 2.5. 



scribed in m ore detail in Szomor u et al.| (|2010) and Szomom 
et al.| ( |2012| ), is different in the sense that these deviations 
are explicitly included in the measurement process. The in- 
trinsic profile is derived by first fitting a Sersic model pro- 
file convolved with the PSF to the observed flux, using a PSF 
constructed from unsaturated bright stars in the image. The 
residuals from this fit are then measured in radially concen- 
tric ellipses which follow the geometry of the best-fit Sersic 
profile. By adding these residuals to the best-fit Sersic profile, 
we effectively perform a first-order correction for deviations 
from the model profile and are able to account for complex 
substructures such as rings and faint outer disks. The result- 
ing profil es closely follow the true int rinsic galaxy profiles , as 
shown in Szomor uet al. | ( [20 10] ) and |Szomoru et al.| ( |2012| ). 

The surface brightness profiles of all galaxies in the high- 
redshift sample are measured in the #435, V606, ^775, Zsso, ^125 
and if 160 filter; this ensures sufficient wavelength coverage to 
accurately measure rest- frame u-g colors. The SDSS galaxy 
profiles are measured in the u, g, r and z bands. Errors in the 
flux profiles are calculated by adding the formal flux errors, 
sky variance, and the estimated error in the sky background 
determination. The radial extent of the profiles is mostly lim- 
ited by uncertainties in the sky background estimation; typ- 
ically, the profiles are accurate out to radii of approximately 
10 kpc. At larger radii the profiles are extrapolated using the 
Sersic model profile. The extrapolated part of the profile typ- 
ically contains ~ 5-10% of the total flux, depending on the 
galaxy profile shape. 

Rest-frame u and g band profiles are derived by inter- 
polating between the observed fluxes at each radius, using 
the SE D-based interpolation package InterRest (Taylor et al. 
2009| ). This package uses a set of template SEDs to interpo 



late between observed fluxes. The resulting rest-frame fluxes 

7 See|http : / /www . mpa-garching . mpg . de/SDSS/DR7/mass_| 
comp . html for a comparison between these masses and masses based on 
spectral indices 
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thus take into account the shape of a galaxy's SED and the 
filter throughputs of both the observed and rest-frame filters. 

3.2. From colors to mass-to-light ratios 

Stellar mass-to-light ratios are estimated using an empir- 
ical relation between W336 - £475 and M*/L g . The u and g 
filters straddle the B aimer and 4000 A breaks; u-g color is 
therefore sensitive to changes in stellar age, dust, and metal- 
licity. Our empirical u-g — M*/L g relation is based on 
FIREWORKS photometry of g alaxies in the CDFS, shown 
in Figure [T] ( Wuyts et al. 



2008 ). In this Figure we plot the 



integrated rest-frame u-g colors and \ogM*/L g of galaxies 
in three redshift bins (grey circl es). Rest - frame colors have 
been calculated using InterRest (Tay lor et al.||2009| ), and the 
mass-to-light ratios are obtained from the SEE) fits described 
in Section |2 The black lines indicate the best- fit linear rela- 
tion at each redshift. The slope of the relation is 1.0 at z ~ 2, 
0.9 at z ~ 1, and 0.7 at z ~ 0. The uncertainty in logM*/L g , 
given a value for u-g, is approximately 0.28 at z ~ 2, and 
0.13 at lower redshifts. 

The reason that these relations exist is due to the fact 
that stellar population variations (i.e., changes in stellar age, 
metallicity, and dust content) all produce r oughly the same e f- 
fects in the u-g — M*/L g plane ( [Bell & de Jong| (|2001|)). 
This is indicated by the arrows in each panel ot Figure IT] 
which all point roughly in the same direction. Due to this de- 
generacy we can not distinguish between dust, age, and metal- 
licity effects, but M*/L g can be determined quite robustly. It 
should also be noted that the axes in Figure [T] are not com- 
pletely independent, since the rest- frame u-g color enters into 
the stellar population fit from which M*/L g is determined. 

3.3. Mass profile derivation 

The process of deriving stellar mass surface density profiles 
is illustrated for four galaxies in Figure [2] The galaxies are 
selected to show a range of color gradients. For each galaxy 
we show, from top to bottom, a rest-frame ubg color image, 
the observed-frame and rest- frame surface brightness profiles, 
the rest- frame u-g color profile, and the resulting stellar mass 
surface density profile. The extent of the PSF half width at 
half-maximum (HWHM) is indicated by grey hatched areas. 

The observed-frame residual-corrected surface brightness 
profiles, shown in grey-scale in the second row, are gener- 
ally of high quality. The profiles deteriorate somewhat in 
the bluer bands for the highest-redshift galaxies. However, 
the profiles that are used for interpolating to rest- frame u and 
g wavelengths (i.e., measured in the bands directly red- and 
blueward of the rest-frame u and g wavelengths) have high 
signal-to-noise. The observed surface brightness profiles, and 
by extension the rest- frame u-g profiles, are generally accu- 
rate out to ~ 10 kpc. 

In the third row of Figure [2] we plot the observed u-g color 
profiles as well as the average u-g colors for the entire galaxy 
(in red and blue, respectively). The color gradients range from 
very steep (A(u - g) / Alogr = -1) to very shallow. At large 
radii (r > 10 kpc) the color profiles become extremely uncer- 
tain. Due to the low surface brightness at these radii, the ratio 
of u band to g band flux is very sensitive to small errors in 
either flux profile. This can result in colors that become un- 
realistically blue or red. We therefore limit the color profiles 
to radii where the error in u-g is smaller than 0.2 dex; the 
u-g color at larger radii is fixed to the value at the threshold 
radius. This should not have a strong effect on the resulting 



stellar mass surface density profile, since the surface densities 
at these large radii are so low that they contribute very little to 
the total mass, even for high M*/L g . 

The resulting mass profiles are shown in the fourth row, 
with the same color coding (blue assuming a constant M*/L, 
red for radially varying M*/L). Effective radii obtained by in- 
tegrating the surface density profiles are shown for both pro- 
files, with the hatched areas indicating the 1 - a errors. The 
difference between assuming a radially constant mass-to-light 
ratio and actually accounting for M*/L variations is clear. Ra- 
dial M*/L variations can lead to half-mass radii that are more 
than a factor 5 smaller than rest-frame g band half-light radii, 
albeit with fairly large errors (first column). Less extreme, 
but still very large, differences of a factor ~ 2 are measured 
with very low errors (second column). These size differences 
will be discussed in more detail in the next Section. The sur- 
face density profiles of all galaxies in our sample are shown 
in Figure [3] and are also given in Table [T] 

4. MASS-WEIGHTED SIZES 

4.1. Structural parameter derivation 

Having constructed stellar mass surface density profiles, 
we focus our attention on galaxy sizes. Half-mass radii are 
amongst the most basic galaxy characteristics and are espe- 
cially important in the context of galaxy evolution, given their 
close connection to the build-up of stellar mass. Half-light 
radii, used in all studies of galaxy size evolution, can be 
strongly influenced by stellar population differences within 
a given galaxy. The surface density profiles described in the 
previous Section provide information on the true mass distri- 
bution of these galaxies. 

Half-mass radii are calculated by integrating the surface 
density profiles out to a radius of ~ 100 kpc. This corresponds 
to ~ 12 arcsec for galaxies at z ~ 1-2, and ~ 90 arcsec for 
the z = galaxies. The errors in the half-mass radii (and other 
parameters derived from the surface density profiles) are es- 
timated in two ways. Firstly, we estimate the uncertainty due 
to errors in the flux profile by varying the u-g color profiles 
within their 1 -a error ranges, deriving mass profiles based on 
these new color profiles, and calculating the resulting range 
of half-mass radii. Secondly, we estimate the effects of PSF 
variations on our size determinations. This is done by rerun- 
ning the entire procedure described in Section [3] using twelve 
different stellar PSFs. The median error in r mass due to these 
sources is 18%. 

In Figure]?] the half-mass radii are plotted against half-light 
radii measured from the rest-frame g band surface brightness 
profiles. The half-mass radii are generally smaller than the 
half-light radii: (log r mass /r g ) w -0.12 ± 0.01. Overall, the 
stellar mass distributions of massive galaxies are more con- 
centrated than their rest-frame optical light distributions by 
about 25%. The overall trend is in agreement with previous 
studies, which have found that galaxies at these redshifts tend 
to show negati ve color gradients, such that their cores are rela- 



tively red (e.g ., van Dokku m et al.|2010"t Szomoru et al.|2011 
|Guo et al.|201ip . 

The median differences between mass-weighted and 
luminosity- weighted parameters are given in Table [2] Mass- 
weighted Sersic indices (obtained by fits to the radial sur- 
face density profiles) are on average ~ 5 - 20% larger than 
luminosity- weighted Sersic indices. It should be noted that 
Sersic indices are much more difficult to accurately constrain 
than effective radii, due to their sensitivity to systematic un- 
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Figure 1. The relationship between rest-frame u-g color and stellar mass-to-light ratio in three redshift bins. Gr ey circles denote in tegrated colors and mass-to- 
light ratios of individual galaxies in the Chandra Deep Field South, based on data from the FIREWORKS catalog ( Wuyts et al. 2008 I. The best-fit linear relations 
are indicated by black lines. The slope of the relation is close to unity at high redshift, and flattens off at lower redshifts. The effects of age, dust extinction and 
metallicity on u-g and M*/L g are indicated by the arrows labeled At, AAy and AZ, respectively. The length of the vectors indicate the shift caused by an 
increase in stellar age from 1 Gyr to 5 Gyr, lAy of dust extinction, and an increase in metallicity from 0.02 to 2.5 times solar metallicity. Increases in stellar age, 
dust content, and metallicity all result in a shift roughly parallel to the best-fit relation, and are thus implicitly included in our empirical relation. 
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Stellar mass surface density profiles. 






ID a redshift rarcsec r kpc 


flu 


logM*/L g 


logS 


(arcsec) (kpc) 


(AB mag arcsec -2 ) (AB mag arcsec -2 ) 


(logM /L ) 


(logM kpc -2 ) 



981 


1 


13 


0.0180 


0.148 


21.470 ±0.0056 


19.442 ±0.0008 


0.671 ±0.1301 


10.552±0.1301 


981 


1 


13 


0.0198 


0.163 


21.509 ±0.0058 


19.499 ±0.0009 


0.655 ±0.1301 


10.514±0.1301 


981 


1 


13 


0.0216 


0.177 


21.547 ±0.0063 


19.552 ±0.0009 


0.642 ±0.1301 


10.479 ±0.1301 


981 


1 


13 


0.0240 


0.197 


21.596 ±0.0083 


19.620 ±0.0011 


0.625 ±0.1302 


10.435 ±0.1302 


981 


1 


13 


0.0264 


0.217 


21.642 ±0.0087 


19.683 ±0.0015 


0.610±0.1302 


10.394 ±0.1302 


981 


1 


13 


0.0288 


0.237 


21.687 ±0.0090 


19.744 ±0.0016 


0.596 ±0.1303 


10.356±0.1303 


981 


1 


13 


0.0318 


0.261 


21.740 ±0.0099 


19.815±0.0018 


0.580 ±0.1303 


10.312±0.1303 


981 


1 


13 


0.0348 


0.286 


21.792 ±0.0124 


19.883 ±0.0019 


0.566 ±0.1305 


10.271 ±0.1305 


981 


1 


13 


0.0384 


0.315 


21.851 ±0.0131 


19.960 ±0.0020 


0.550 ±0.1305 


10.224 ±0.1305 


981 


1 


13 


0.0426 


0.350 


21.919±0.0139 


20.044 ±0.0023 


0.535 ±0.1306 


10.176±0.1306 



Note. — This Table is published in its entirety in the electronic edition of this Paper, and can also be downloaded from http: 
|//www . strw . leidenuniv . nl/~szomoru/| A portion is shown here for guidance regarding its form and content. 
a FIREWORKS ID jWuyts et al.|20Q8) for 0.5 < z < 2.5 sources, NYU-VAGC ID (Blanton et alfe OQSi for z < 0.5 sources. 



certainties; errors on our mass- weighted Sersic indices are ap- 
proximately twice as large as errors on our mass- weighted ef- 
fective radii. 

We also give the difference between total masses calculated 
by integrating the stellar mass surface density profiles and to- 
tal masses from integrated photometry. We find that, on aver- 
age, masses from resolved photometry are higher than those 
based on integrated light: the median difference is 5 - 10% 
at high redshift, and ~ 30% at z = 0. Our low-redshift result 



is in agreement with results presented by Zibet ti et al.| ([2009), 
who derive resolved stellar mass maps of nine nearby galaxies 
with a range of morphologies. These authors find that mass 
estimates from integrated photometry may miss up to 40% 
of the total stellar mass compared to estimates obtained by 
summing resolved mass maps, due to dusty regions being un- 
derrepresented in the total flux of galaxies. 

This has consequences for integrated galaxy colors, which 
are often used as a proxy for star formation activity. Dusty, 
red regions with low flux but high mass densities will be un- 
derrepresented in colors based on integrated galaxy photome- 
try. We estimate this effect using total rest- frame u-g colors. 



Mass- weighted colors are calculated as follows: 



(u-g) n 



fZ(r)dr ' 



(1) 



where f u (r), f g (r) and E(r) are the radial u and g band flux 
density and radial stellar mass surface density profiles, respec- 
tively. On average, mass-weighted u-g colors are redder than 
luminosity- weighted u-g colors by 0.04-0.06 magnitudes, 
indicating that the SSFRs implied by luminosity- weighted 
u-g colors are slightly overestimated compared to those im- 
plied by mass- weighted colors. 

4.2. Sources of uncertainty 

The stellar mass surface density profiles are generally not 
very sensitive to the extrapolation of the color profiles towards 
larger radii, due to the fact that the total flux at r > 10 kpc is 
very low. Furthermore, this source of uncertainty is explicitly 
included in our errors, since they are estimated by varying the 
u-g color that is used for the extrapolation. Uncertainty in 
the shape of the color profile within the PSF HWHM is more 
difficult to quantify, and may be especially important for the 
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Figure 2. Illustration of the conversion of observed surface brightness profiles to stellar mass surface density profiles, for four sample galaxies at 0.5 < z < 2.5. 
Top row: color images, composed of rest-frame W336, #433 and #475 images. Second row: observed-frame residual-corrected surface brightness profiles (shown 
in greyscale), with rest-frame u and g band profiles overlaid in purple and yellow, respectively. Third row: measured u-g color profiles and average u-g colors 
(in red and blue, respectively). The grey hatched area indicates the PSF HWHM. Bottom row: resulting surface density profiles, with color coding corresponding 
to the color profiles in the third row. Effective radii are shown for both profiles; the red and blue hatched areas indicate the 1 -a errors. Observed- and rest-frame 
photometry is generally of very high quality, and provides accurate u-g color profiles out to ~ 10 kpc. The colors at the largest radii, where the flux profiles are 
not reliable, are extrapolated from the last well-measured color. Half-mass radii are in some cases more than 50% smaller than half-light radii. 



Table 2 

Differences between mass-weighted and luminosity-weighted parameters. 







log r mass - 


l0g?g 


log n mass - 


-k)g% 


log Mmass- 


logM/ 


(u g)mass 


(u-gf g 


median 


ONMAD 


median 


ONMAD 


median 


ONMAD 


median 


ONMAD 


z 


= 


-0.12 ±0.01 


0.08 


0.06 ±0.03 


0.31 


0.12±0.01 


0.14 


0.059 ±0.004 


0.048 


0.5 < 


z< 1.5 


-0.14±0.01 


0.11 


0.09 ±0.02 


0.20 


0.04 ±0.01 


0.06 


0.044 ±0.006 


0.048 


1.5 < 


z<2.5 


-0.10 ±0.02 


0.13 


0.02 ±0.03 


0.19 


0.02 ±0.01 


0.07 


0.043 ±0.007 


0.051 



Note: the mass subscript indicates parameters derived using the true mass profiles (i.e., with radially varying M*/L), while the g subscript indicates 
parameters derived using profiles that assume a constant M*/L profile (i.e., equivalent to the rest-frame g band profiles). 

a logM mass is the total stellar mass derived by summing the resolved mass density information. logMg is the total stellar mass based on integrated 
photometry. 

b {u-g) mass is the mass-weighted color (see Equation|4.l|, and {u-g) g is the luminosity- weighted color. 
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z = 0.5<z<1.5 1.5<z<2.5 
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Figure 3. Stellar mass surface density profiles of galaxies with M* > 1O 1O 7 M0 at z = 0, z ~ 1, and z ~ 2. Individual profiles are shown in blue and red for 
starforming and quiescent galaxies, respectively. The approximate PSF HWHM in each redshift bin is indicated by the hatched region. The profiles are plotted 
up to the radius where the errors become significant; they are accurate out to ~ 10 kpc, and down to surface densities of around 1O 7 5 M0 kpc -2 . 




Figure 4. Half-mass radii versus rest-frame g band half-light radii, for galaxies with M* > 1O 1O 7 M0. The approximate PSF HWHM in each redshift bin 
is indicated by the hatched region. At all redshifts the half-mass radii of galaxies are on average 25% smaller than their half-light radii. Values as low as 
log r ma ss/r g = -0.7 are found. 
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Figure 5. Comparison of half-mass radii derived using CANDELS GOODS- 
South data and HUDF data, for the galaxies with 0.5 < z < 2.5 and M* > 
10 10 - 7 that are found in both datasets. The two datasets agree within the 
errors. 

smallest galaxies (r g < 1 kpc). In order to estimate the stabil- 
ity of our derived sizes we have explored several alternative 
approaches: one in which the inner u-g color is kept fixed to 
the integrated value within 1 kpc; and one in which the best- 
fit linear color gradient is extrapolated from the PSF HWHM 
inward. Neither approach alters our results. 

The measured shape and radial extent of a galaxy's sur- 
face brightness profile are sensitive to the depth to which 
the galaxy is imaged. A lack of imaging depth can result 
in errors in sky background estimation, as well as a portion 
of the galaxy's emission being lost in the background noise. 
Generally, low-S/N data is likely to introduce systematic ef- 
fects, such that measured sizes and Sersic indices are smaller 



tha n the true values (e.g., |Trujillo et al. 2006, Williams et 



lar][2010). As a consistency check we therefore compare 
our results to those based o n ultradeep optical and NIR data 
acquired ov er the HUDF (Beckwith et al.||2006[ |Bouwens 
et al. |2011| ). These data are ~ 2 magnitudes deeper in the 
than thi 



NIR than the CANDELS data and should therefore be unaf- 
fected by these surface brightness effects. We find 1 1 galax- 
ies with M* > 1O 1O - 7 M and 0.5 < z < 2.5 which are imaged 
by both CANDELS and HUDF. The HUDF WFC3 data of 
one of these galaxies shows some artificial background fea- 
tures, which are the result of the background subtraction pro- 
cess; this galaxy is therefore excluded from the comparison. 
The HUDF imaging is processed in exactly the same way as 
described in Section |3] for the CANDELS imaging. The re- 
sulting half-mass radii are compared to the CANDELS half- 
mass radii in Figure [5] The correspondence is good; HUDF 
and CANDELS measurements of r mass lie within la of each 
other, and we find no systematic offset. Thus, based on this 
subsample of galaxies we conclude that the CANDELS data 
are sufficiently deep for our purposes. 

4.3. Evolution with redshift 

Figure [4] indicates that mass distributions are on average 
more concentrated than light distributions at all redshifts. We 
now investigate the redshift evolution of this concentration 
difference in more detail. In Figure [6] we plot r mass /r g as a 
function of redshift for all galaxies in our sample, as well as 
for quiescent and starforming galaxies separately. Quiescent 
galaxies are defined to have SSFR < 0.3/r#, where tu is the 



Hubble time. Individual galaxies are indicated by small cir- 
cles, and median values for each redshift bin are shown as 
large, darker circles. The error bars on the median values in- 
dicate the scatter of the distribution. 

Overall, the difference between mass-weighted and 
luminosity- weighted radius does not seem to evolve with red- 
shift: log r mass /r g = -0.12 ± 0.01 at z = 0, -0.14 ± 0.01 at 
z ~ 1, and -0.10 ± 0.02 at z ~ 2 (see also Table J2j. Simi- 
larly, the values for starforming galaxies are consistent with 
no evolution with redshift. For quiescent galaxies, there is a 
hint of decreasing values of log r mass /r g at higher redshifts, al- 
though the large scatter at z ~ 2 means that these results are 
also consistent with zero redshift evolution. 

The lack of evolutio n with redshift agrees w ith the color 
gradients presented in Szomoru et al.| (2011|). This study 



showed that the radial color gradients in galaxies with 10 iU < 
M*/M < 10 11 are nearly constant between < z < 2.5. This 
seems counterintuitive, since one would expect bulge growth 
in galaxies to result in steeper color gradients at lower red- 
shift. However, it should be stressed that the comparison in 
Figure |6]is between galaxies of the same mass. It is therefore 
not a comparison between high-redshift galaxies and their 
descendants, but rather a comparison between high-redshift 
galaxies and their low-redshift analogs. Thus, z ~ 2 galax- 
ies have similar color gradients as low redshift galaxies of the 
same mass. This holds for starforming and quiescent galaxies 
seperately, and for the entire population as a whole. 



4.4. Dependence on galaxy parameters 
The transition of galaxies from the starforming to the pas- 



of their structure and morphology (e.g., Kauffmann et al. 


2003 Toft et al. 


2007; |Franx et al.||2008[ |Bell 20081 Ivan 


Dokkum et al. 


2011 ISzomoru et al. 201 1| Wuyts et al. 


2011; Bell et al 


. 2012). The growth of a bulge, in partic- 



are quenched. We therefore expect some degree of correla- 
tion between r mass /r g and the galaxy parameters that correlate 
with star forming activity, such as SSFR, size, Sersic index, 
and effective surface density. We investigate these correla- 
tions in Figure [7] Individual galaxies are indicated by small 
light-colored circles, and the running medians in each panel 
are indicated by darker lines. 

Somewhat surprisingly, the median value of r mass /r g is 
close to constant as a function of galaxy parameters in all 
redshift bins. There is some evidence for a trend between 
r mass /r g and galaxy morphology, such that r mass /r g is smaller 
for galaxies with low SSFR, small sizes and high Sersic in- 
dices. However, this trend is very weak. Overall, the cor- 
relation coefficients at z ~ 1 and z ~ 2 are very low and not 
significant (r ~ 0.01, p ~ 0.4). The exceptions to this are the 
correlations between \ogr mass /r g and log SSFR and n serS i C , at 
z ~ 2; although these are also quite low, they are significant 
(r ~ 0.35, p < 0.01). Similarly, the correlation coefficients 
at z = are low, but significant (r ~ 0.2-0.3, p < 0.01). We 
conclude that the difference between half-mass size and half- 
light size correlates very weakly with galaxy structure and star 
forming activity. 

Although the average difference between half-mass and 
half-light radius is remarkably constant, for some galaxies 
the difference is up to a factor eight. These large differences 
warrant closer inspection. In Figure [8] we show images and 
profiles of a number of galaxies witn r mass /r g < 0.5. Four 
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All galaxies Quiescent Starforming 




redshift redshift redshift 



Figure 6. The ratio of half-mass size to half-light size, r mass /r g , as a function of redshift, for galaxies with M* > 10 107 . From left to right are plotted all galaxies, 
quiescent galaxies, and starforming galaxies. The large symbols indicate the running median, with the error bars indicating the 1-cr scatter of the distribution. 
Our results indicate that, at fixed mass, galaxies have similar r mass / r g values at all redshifts between z = and z = 2.5. 




Figure 7. The correlations between r maS s/r g and, from left to right: total stellar mass, specific star formation rate, and structural parameters measured in the 
rest-frame g band (effective surface density, effective radius, and Sersic index). Individual galaxy measurements are indicated by circles, and running medians 
are indicated by the thick lines. The error bars on the running medians indicate the 1-cr scatter. Each row corresponds to a redshift bin. r mass /r g correlates 
weakly, if at all, with starforming activity at both high and low redshift. 



representative galaxies have been selected from each redshift 
bin, and are plotted in order of decreasing SSFR. The flux 
profile errors of the galaxies shown are typical for the galax- 
ies in this sample. Nearly all of these galaxies are dominated 
by large central bulges, but also contain extended, blue disks. 
These large galaxies have very well-measured surface bright- 
ness profiles; errors are low, and the profiles are measured 
accurately out to very large radii. We are therefore confident 
that the color gradients, although very steep, are real. 

This group of galaxies demonstrates the importance of ob- 
taining resolved stellar mass information. Measurements 
made at rest-frame optical wavelengths severely underesti- 
mate the relative importance of the central bulge component 
in these galaxies, due to the high luminosity of the blue disks 



surrounding them. As a result, for a sizeable fraction (~ 10%) 
of massive galaxies the structure as measured at rest-frame 
optical wavelengths severely misrepresents the mass distribu- 
tion in these galaxies. 

This population of large, bulge-dominated disk galaxies is 
very interesting. They have stellar masses up to ~ 1O 11 M 
and are distributed across the entire range of measured SS- 
FRs. Moving from high to low SSFR, we can see a progres- 
sion from Sa or Sb-like systems to objects that look more like 
SOs and elliptical galaxies. This is most evident at low red- 
shift, where galaxies have large angular sizes. The color dif- 
ference between the central bulge and extended disk gradually 
decreases from high- to low-SSFR galaxies. In other words, 
it seems that we are observing these galaxies during the phase 






1 10 1 10 1 10 1 10 

radius (kpc) radius (kpc) radius (kpc) radius (kpc) 



Figure 8. A representative selection of galaxies with r mass /r g < 0.5, ordered by SSFR. All galaxies are extended, disk-like galaxies with large central bulges. 
Measurements in the rest-frame g band result in severely overestimated sizes compared to the half-mass sizes. The color gradients are measured with high 
precision; the large differences between r mass and r g are therefore not the result of measurement errors. 
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where a central bulge has recently formed and star formation 
is turning off. It is unclear why other galaxies with similar SS- 
FRs do not contain such large bulges; this may depend on the 
dark matter halos in which they are located, or may perhaps 
be due to different quenching processes. 

5. SUMMARY AND CONCLUSIONS 

We have presented the first consistently measured stellar 
mass surface density profiles for individual galaxies at red- 
shifts between z = and z = 2.5. These profiles have been de- 
rived using an empirical relation between rest- frame color and 
stellar mass-to-light ratio. This simple method does not yield 
the same detailed information regarding, e.g., stellar ages 
and dust content as resolved SED-fitting techniques. How- 
ever, it is robust to variations in stellar population properties; 
changes in stellar age, metallicity, or dust content shift galax- 
ies roughly along the empirical relation, and are thus implic- 
itly included in our M*/L determination. 

The key advantage of this study compared to previous work 
at hig h redshift (e.g., |Wuyts et al.|2012||Lany on-Foster et al 



2012) is the fact that the profiles presented here are decon 



volved for PSF smoothing. This is crucial for measurements 
of high-redshift galaxy structure, since these distant and of- 
ten physically small galaxies have angular sizes t hat are in 
many cases comparable to the HST PSF size (e.g., [Dad di et 
al. 120051 jTnljillo et al.|2006[IToft et al.|2007[rv an Dokkumet^ 
al.|2008[ [Szomoru et al.|2010l|Cassata et al.|2010n Szomoru et 



^11 

al. 2012). Since our surface density profiles are derived from 
deconvolved surface brightness profiles, they can be used to 
correctly measure structural parameters such as sizes and Ser- 
sic indices. 

The considerable depth of the data used in this study al- 
lows us to probe galaxy fluxes and colors out to large radii. 
The robustness of the resulting structural parameters has been 
tested using ultradeep data taken over the HUDF, which over- 
laps with the CANDELS GOODS-South field. A compari- 
son of galaxy parameters derived using the two datasets con- 
firms that our results are not systematically affected by surface 
brightness effects. 

We have shown that the half-mass radii of galaxies between 
z = and z = 2.5 are on average 25% smaller than their rest- 
frame optical half-light radii. There is significant scatter in 
this size difference, with some galaxies having half-mass radii 
that are almost an order of magnitude smaller than their half- 
light radii. We find that, on average, this size difference does 
not vary with redshift for galaxies at fixed mass. This holds 
for the population as a whole, as well as for the quiescent and 
starforming subpopulations seperately. This is an interesting 
result, as it implies that z ~ 2 galaxies have similar color gra- 
dients as their low-redshift analogs, despite the fact that these 
low-redshift galaxies formed at a different epoch, and perhaps 
through very different formation mechanisms. 

There does not seem to be a strong correlation between 
galaxy morphology or star forming activity and the difference 
between half-mass size and half-light size. However, we do 
find that the galaxies with the most extreme size differences 
are almost all extended disk galaxies with very prominent cen- 
tral bulges. These galaxies range from strongly starforming to 
almost completely quiescent, and may represent a short tran- 
sitional phase during which the central bulge is prominent and 
the starforming disk is very young. 

We note that there is significant scatter around the empiri- 
cal relation used to convert u-g color to M*/L. Although age, 
metallicity, and dust effects are implicitly included in this re- 



lation, their effects are not exactly parallel to the relation. This 
could lead to a systematic underestimate of the mass-to-light 
ratios in galaxy regions that are relatively metal-rich or old, 
and an overestimate in regions that are relatively metal-poor 
or young. Since the central regions of galaxies generally con- 
tain older stars, the inferred M*/L gradients would then be too 
shallow. This effect is likely small, since the vectors shown in 
Figure |T] are close to parallel, and the magnitude of the shift 
is smallfor moderate stellar population differences. 

A possible concern is regions with very high dust content, 
which could potentially obscure stellar light to such a degree 
that our u-g — M*/L conversion becomes ineffective, sim- 
ply because all stellar light is extincted. Very high central dust 
concentrations may result in an underestimate of the inner 
mass content of galaxies, and an overestimate of their half- 
mass radii. Such effects are a significant source of system- 
atic uncertainty in our analysis, and can only be quantified by 
measuring the light reemitted by dust. This requires infrared 
imaging at HST resolution or better, which will become pos- 
sible in the near future using ALMA. 

Inside-out g a laxy growth, as described by, e.g., |van| 
Dokkum et al. ( 2010| ), implies that the growth of the most 
massive galaxies since z ~ 2 is largely due to material be- 
ing accreted onto the outer regions of these galaxies. The 
cores of massive galaxies likely formed in short, violent bursts 
at higher redshift, and should therefore have star formation 
histories and stellar populations that are quite different from 
those in the outer regions. The results presented in this Pa- 
per broadly agree with such a picture; the central regions of 
massive galaxies are redder, and therefore likely older, than 
the outer regions. Using the method presented in this Pa- 
per we cannot, however, disentangle dust, age and metallic- 
ity gradients; nor can we constrain the star formation histo- 
ries within our galaxies. First steps towards a better under- 
standing of stellar population variations within high-redshift 
galaxies have been made by several author s. Results based 
on photometry of ear ly-type galaxies (e.g., |Guo et aL]|2011 



Garg iulo et al. 12012] ) indicate that stellar age and metallic 



ity are the dominant drivers of radial color gradients. Stud 
ies based on spectroscopic measure ments of gravitationally 
lensed high-redshift galaxies ( e.g., |Cresci et al."||2010t Jones 
etaLl[20i0l |Yuan et al ][20TT] [Queyrel et al.||2012| | J ones el 
al.||2012| ) nave shown that most of these galaxies have nega 



tive metallicity gradients. These results seem to be roughly 
consistent with each other, but are based on very differently 
selected, and rather small, galaxy samples. A broader, more 
in-depth analysis of radial stellar population variations, for a 
well-defined sample of starforming galaxies as well as qui- 
escent galaxies, could provide valuable insights into the pro- 
cesses which have shaped the structure of galaxies today. 
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